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Summary 


Slemrod’s  reseaxch  in  1990-91  centered  on  two  issues:  (1)  the  kinetics  of  co.ig  nation 
processes,  (2)  behavior  of  discrete  velocity  models  in  the  kinetic  theory  of  gases.  In  the  first 
area  Slemrod  hcis  (a)  given  a  new  method  for  solving  the  special  class  of  coagulation  equa¬ 
tions  which  exhibit  gelation  and  (b)  derived  and  proved  existence  of  similarity  solutions 
for  coagulation  equations  with  diffusion.  In  the  second  area  Slemrod  has  used  his  “relaxed 
invariance  principle”  method  to  prove  weak  decay  to  equilibrium  for  the  Broadwell  model 
of  gas  dynamics  in  the  case  of  specularly  reflective  boundary  conditions. 
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In  the  last  year  Slemrod's  reseaxch  has  centered  on  two  issues:  (1)  the  kinetics  of 
coagulation  processes,  (2)  behavior  of  discrete  velocity  models  in  kinetic  theory  of  gases. 
Below  an  outline  of  progress  in  each  area  is  given. 


1.  Kinetics  of  Coagulation  Processes 

Models  of  cluster  growth  appear  in  a  wide  variety  of  applications.  One  well  known 
exaunple  is  the  Smoluchowski-Flory-Stockmayer  theory  of  gelation  where  it  was  found  that 
all  concentrations  decrease  in  time  and  the  quantity  representing  mass  density  is  conserved 
for  only  a  finite  time  cifter  which  it  decreases. 

The  models  themselves  are  coupled  infinite  systems  of  ordinary  differential  eauations: 

Cj{t)  >  0,  j  =  1,2, ...  denote  the  expected  numbers  of  clusters  consisting  of  j  particles 
per  unit  volume.  The  discrete  coagidation-fragmentation  equations  are 

^  h  ^ 

(1) 

^  ^j,k^j+k]  Ibr  j  =  1,2... . 

The  coagulation  rate  aj^k  and  fragmentation  rate  bj^k  are  non-negative  constants  with 

^j,k  ~  ^k,ji  ^j,k  ~  ^k,j' 

In  paper  [1],  Slemrod  considered  the  case  of  pure  coagulation  for  which  byk  =  0  and  the 
special  case  ajk  =  jk.  In  this  case  the  Smoluchowski-Flory-Stockmayer  system  becomes 

J  1  °° 

~  2  -  jcj  kck  .  (2) 

^  k=l  k=l 


For  simplicity  monodisperse  initiaJ  data  was  chosen  Ci(0)  =  1,  Cj(0)  =  0,  2  <  j  <  oo. 

A  solution  of  (2)  has  been  given  by  McLeod  (J.  B.  McLeod,  On  an  infinite  set  of 
non-linear  differential  equations.  Quarterly  J.  Math.  Oxford  Ser  (2),  13  (1962),  119-128). 
McLeod’s  solution  was  valid  on  0  <  t  <  1.  It  was  based  on  the  fact  that  on  0  <  t  <  1  the 

OO 

mass  density  p{t)  =  ^^kck{t)  is  conserved  (and  equals  one).  However  as  McLeod  noted 

fc=i 

the  desired  conservation  of  density  breaks  down  for  t  >  1  eind  his  solution  is  not  valid  past 
the  critical  “gelation”  time,  i  =  1. 

The  first  resolution  of  the  problem  was  provided  by  Leyvraz  and  Tschudi  (F.  Leyvraz 
and  H.  R.  Tschudi,  Singularities  in  the  kinetics  of  coagulation  processes,  J.  Phys.  A.  14 
(1981),  3389-3405).  Their  method  was  to  introduce  the  generating  function  G{z,t)  = 

OO 

where  =  jCj(t)exp(j  f*  p(z)dz).  A  computation  shows  G  satisfies  the 

k=l 

qucisi-lineau-  hyperbolic  partial  differential  equation 


dG  .dG 
dt  ~^^dz 
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0  <  2  <  1,  t  >  0 


(3) 


with  initial  data  G(2,0)  =  z.  They  integrated  (3)  via  the  method  of  characteristics  and 
then  recovered  Cj{t).  The  end  result  (after  long  computations)  was  that 


.  .  ^  1 

■  (i-1)! 


-it 


t  >  1 

0  <  /  <  1 


In  this  case  the  density  p(t)  satisfies 


(4) 


p(t)  =  1  0  <  t  <  1 

=  —  1  <  <  <  oo 

i 

indicating  a  decrease  in  density  after  the  “gelation”  time  t  =  \. 

In  [1]  Slemrod  reconsidered  resolution  (2)  based  on  McLeod’s  origincd  ideas  with  some 
subtle  changes.  He  was  able  to  show  the  density  p{i)  satisfies  the  singular  ordinary  differ¬ 
ential  equation 

which  immediately  yields  (5)  and  thus  (4).  Hence  no  recourse  to  complicated  use  of 
characteristics  and  generating  functions  is  needed  and  a  clear  picture  of  the  evolution  of 
the  density  is  provided. 

In  a  second  paper  Slemrod  [2]  following  ideas  of  Binder  (K.  Binder,  Theory  for  the 
dynamics  of  clusters,  II.  Critical  diffusion  in  binary  systems  and  the  kinetics  of  phase 
separation,  Phys.  Rev.  B  15  (1977),  4425-4447),  considered  the  incorporation  of  diffusion 
in  the  classical  discrete  coagulation-fragmentation  system  (1).  Again  the  restriction  was 
made  to  pure  coagulation  {bjk  =  0)  but  now  the  added  difficulty  of  diffusion  (ais  given  by 
Stokes’  law)  is  considered. 

The  equations  now  are  for  qucintities  cy(i,  f)  >  0,  j  =  1,2...  which  denote  the  expected 
number  of  clusters  of  j  pzirticles/unit  volume  at  position  i  e  IR"(n  =  1, 2, 3)  at  time  t  >  0. 
The  discrete  coagulation-diffusion  equations  are 


dt 


DjAcj 


j  J-l  oo 

9  ^j-k,kCj-k<^k  ~  ftj.kCjCk  j 
*=1  k=l 

i  =  l,2,...; 


(7) 


Dj  ~  j~^l^  (Stokes’  law). 

In  [2]  Slemrod  considered  a  variety  of  problems:  Existence  of  steady  states,  finite  time 
gelation,  amd  diffuse  interfaces.  A  particularly  interesting  point  is  that  system  (7)  admits 
an  exact  similairity  form  Cj{R,t)  =  where  ^  =  R/yt,  R  =  lx|.  The  functions  4>j{0 

satisfy  the  infinite  systems  of  ordinary  differential  equations 


-I — ^ — 4>'j)  +  +  <i>j  — 

^  7-1  oo 

=  -9  ^a]-k,k<i>j-k<f>k  +  Cj^aj^k<f>j<i>k  ■ 
k=i  k=i 


(8) 
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In  [2]  Slemrod  also  provided  some  preliminary  existence  theorems  for  solvability  of  (8). 


2.  Discrete  Velocity  Models  of  Gases 

Discrete  velocity  models  are  finite  velocity  approximations  of  the  Boltzmann  equa¬ 
tion.  They  provide  a  mathematicad  simplification  for  the  true  Boltzmann  equation  hence 
often  allowing  more  complete  analysis  than  can  presently  be  done  for  the  true  Boltzmann 
equation.  One  such  discrete  velocity  model  is  the  Broadwell  model  which  describes  a  gas 
of  pairticles  with  identical  masses  moving  along  three  perpendicular  coordinates  with  same 
speed  c.  Results  of  a  particular  collision  have  the  same  probability  with  only  binary  colli¬ 
sions  allowed.  Let  Ni  =  Ni{x,y,  z^t),  i  =  1,2, ...,  6,  denote  the  density  of  particles  moving 
in  the  six  allowed  directions.  Then  iVi,...,jV6  satisfy  the  Broadwell  system 


dNi 

A 

dNi 

dt 

*r 

""  dx 

dN2 

dN2 

dt 

dx 

dN2 

dN2 

dt 

+ 

"  dy 

dN4 

dNi 

dt 

"  dy 

dNs 

t 

dNi 

dt 

i" 

""  dz 

dNe 

dNi 

dt 

dz 

=  (TiNiN2  +  NsN6-2N3Ni), 


==<r{NiN2  +  N3Ni-2NM, 


where  a /2c  is  the  cross  section  for  binairy  collisions. 

For  flows  which  are  independent  of  y,  z  and  for  which  Nj  =  =  N5  =  the  above 

six  velocity  model  reduces  to  the  simpler  form 


^-c^=2c(Ni-N,N,),  (10) 

^  =.7(W.JV,-JV|). 

An  open  problem  concerning  (10)  has  been  the  asymptotic  behavior  (as  t  — >  00)  of 
(10)  with  initial  conditions 

iVi(i,0)  =  iVio(x),  iV2(r,0)  =  N2o{x),  N^ix^O)  =  N,o{x),  0  <  i  <  1  (11) 

and  boundzuy  conditions  of  specular  reflection  type 

Ni{x,t)  =  N2[x,t)  for  t  >  0,  1  =  0,1.  (12) 
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In  [3]  Slemrod  gave  a  resolution  fo  this  problem  via  his  “relajced”  invariance  principle  based 
on  dynamical  systems  in  infinite  dimensional  state  spaces  and  Young  measures.  The  main 
result  is  that  as  /  — oo  the  unique  solution  of  (10),  (11),  (12)  tends  to  traveling  waves 
/i(x  —  cf),  /2(x  +  ct),  /3(x)  satisfying  the  collision  inequality 

/i(x  -  cf)/2(x  +  cf)  < /3(x)^  a.e.  xe[0, 1],  f>0. 

Convergence  to  /i(x  —  ct),  /2(x  +  ct),  fzix)  is  in  the  weaJc*  topology  of  an  appropriate 
L\logL\  Banach  space. 
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